A large age for the pulsar B 175 7— 24 from an upper limit on 
its proper motion 
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The "characteristic age" 1 of a pulsar usually is considered to approximate its true age, but this assump- 
tion has led to some puzzling results, including the fact that many pulsars with small characteristic 
ages have no associated supernova remnants. 2,3 The pulsar B1757 24 is located just beyond the edge 
of a supernova remnant 4 6 ; the properties of the system indicate that the pulsar was born at the centre 
of the remnant, but that it has subsequently overtaken the expanding blast-wave. 5 8 With a character- 
istic age of 16,000 yr, this implies an expected 8 proper motion by the pulsar of 63 80 milliarcseconds 
(mas) per year. Here we report observations of the nebula surrounding the pulsar which limit its 
proper motion to less than 25 mas yr -1 , implying a minimum age of 39,000 yr. A more detailed 
analysis argues for a true age as great as 170,000 yr, significantly larger than the characteristic age. 
From this result and other discrepancies associated with pulsars, we conclude that characteristic ages 
seriously underestimate the true ages of pulsars. 

The radio supernova remnant G5.4-1.2 is shown in Fig 1: it has an approximately circular radio morphology 
of diameter ~ 35', but is much brighter and has a significantly flatter spectral index along its western edge. 6 ' 8-10 
Approximately 5' beyond this bright western rim is the 125-ms pulsar B1757-24, 4-6 whose rotational frequency, v, 
and frequency derivative, v, can be used to derive a characteristic age, r = z//2z> rj 16 kyr. Provided that a pulsar's 
initial frequency is much larger than v, that the pulsar slows down by magnetic dipole braking, and that its magnetic 
dipole moment and moment of inertia do not evolve with time, then r m t p , where t p is the pulsar's true age. 1 

Surrounding the pulsar is the radio nebula G5. 27-0. 90, whose properties are consistent with it being powered by 
the pulsar's relativistic wind. 6,8 G5. 27-0. 90 has a cometary morphology with a tail pointing back towards the centre 
of G5.4-1.2, an appearance which strongly suggests that the pulsar is travelling away from the centre of G5.4-1.2 
at high velocity. The characteristic age of B1757-24 and its positional offset from the centre of G5.4-1.2 together 
imply a proper motion of 63-80 milliarcsec yr -1 (ref. 8). For a distance to the system of 5 kpc, 8 this corresponds to 
a transverse velocity for the pulsar in the range 1500-1900 km s , significantly higher than typical pulsar velocities 
of -500 km s- 1 (rcfs 11, 12). 

To look for this expected proper motion, we have undertaken multi-epoch observations of G5. 27-0. 90 with the Very 
Large Array (VLA). Our observations are summarised in Table 1; the two epochs are separated by 6.7 yr. The two 
data-sets were processed in identical fashion, with results shown in Fig 2. If G5. 27-0. 90 is associated with the pulsar, 
then the expected shift of the nebula between these two epochs is 422-526 milliarcsec in a westwardly direction; it is 
immediately apparent, even by eye, that motion of this magnitude is not observed. The two epochs can be compared 
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in a more quantitative fashion using the approach described in the legend to Fig 2; we find an offset between the two 
epochs (in the sense [Epoch 2]-[Epoch 1]) of A(RA) = +39 ± 41 milliarcsec and A(Dec) = — 9 ± 34 milliarcsec, where 
contributions from both the noise in the images and estimated errors in the phase calibration of the data have been 
combined in quadrature. We can thus rule out the proper motion which would be expected if t p = t at the 11-14 u 
level, independent of the distance adopted to either the pulsar or the SNR. 

We now consider various explanations for this surprising discrepancy, and show that the most reasonable explana- 
tion is that t p 3> t. In the following discussion, we only consider proper motion in a westward direction (as implied 
by the morphology of G5. 27-0. 90), and adopt a 5-cr upper limit on any shift between epochs of 166 milliarcsec, corre- 
sponding to a maximum proper motion of 24.8 milliarcsec yr -1 . For a distance d — 5d$ kpc and a projected velocity 
of the nebula V n — 1000w„ km s" 1 , this lets us derive an upper limit v n < 0.59^5. 

One possible explanation for this unexpectedly low velocity is that the pulsar is moving at V p <~ 1500 km s _1 w 
2.5v n , but its associated nebula, sensitive to its surroundings, does not track the ballistic motion of the pulsar on 
a one-to-one basis. The distance, r w , between the pulsar and the head of the bow-shock is determined by pressure 
balance between the pulsar wind and the ram pressure produced by the pulsar's motion, such that (ref. 6, 13): 

vE/^rl c = ^pV*, (1) 

where p is the mass density of the ambient medium and 77 is the fraction of the pulsar's spin-down luminosity, 
E = —A-k 2 Ivv = 2.6 x 10 36 erg s" 1 , which goes into powering its relativistic wind (77 < 1). Using the pulsar's position 
as measured in late 1999 (see legend to Fig 2), we find that r w = 0.036^5 pc for the second epoch. Assuming that 
v p = 2.5v n (where V p — 1000u p km s" 1 ), we can then infer that r w must have decreased by at least 20% between 
Epochs 1 and 2. A decrease in r\E by a similar factor can produce this effect. However, comparison of pulsar timing 
parameters from the time of Epoch 2 (A. G. Lyne, private communication) with those published before Epoch 1 (ref. 
5) shows no significant change in E, while a change in r\ seems unlikely given that changes in the efficiency of pulsar 
winds, if they occur at all, most likely occur on time-scales of <~ 10 4 — 10 6 yr. 14,15 

A reasonable alternative is that the pulsar is propagating into a density gradient so that p is increasing as a function 
of time, causing a decrease in r w . However, an increasing density should result in a change in the morphology or 
surface brightness of the leading edge of G5. 27-0. 90. In fact, the structure and brightness of the head of G5. 27-0. 90 
is remarkably similar between the two epochs, and also in archival 6 cm VLA data taken at various epochs between 
1985 and 1999. There is thus no evidence from the morphology of G5. 27-0. 90 for a change in ambient conditions 
between the two epochs. 

It thus seems most likely that r w has not changed between the two epochs, implying that V p = V n -C 1500 km s _1 . 
In this case, demanding that t p ~ t would mean that the pulsar had only travelled 7' over its lifetime, inconsistent 
with its separation of 16.1' — 20.6' from the apparent centre of the SNR 8 and thus arguing against a physical associ- 
ation between the two objects. While such chance coincidence along the line-of-sight is not impossible, 16 the strong 
morphological evidence for an association 5-8,10 leads us to conclude that t p > 39 — 50 kyr » r. This result depends 
only on assuming that V p = V n , and that the pulsar was born near the SNR's apparent centre. 

Another approach is to let the age of the system, the ambient density and the transverse velocity of the pulsar all 
be free parameters. We can then solve for these unknowns using the three characteristic length scales of the system, 
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namely the radius of the SNR (R s = 25^5 pc), the stand-off distance of the bow-shock (r w = 0.036e? 5 pc) and the 
offset of the pulsar from the SNR's centre (R p = [23.4 — 30.0]d5 pc). To determine how the unknown quantities 
relate to R s , we first need to determine the SNR's evolutionary state. Using our upper limit v p < 0.59d 5 and the 
value measured for r w , Equation (1) implies a particle density n\ > n = 0.23r/<if cm -3 . Assuming that the SNR 
is expanding into a similar environment, it will be in the radiative (or "snowplough" ) phase of its evolution 7 ' 17 if 
t p > tji = 13i?5i 14 rj^ 4 / 7 kyr, where E 5 il0 51 erg is the energy in the supernova explosion and assuming solar ISM 
abundances. The lower limits we have already computed on m and t p ensure that the SNR is indeed in the radiative 
phase; the radius of the SNR is then given by: 7 ' 17 

R s = 7.0E°? 21 ni°- 257 4- 3 pc, (2) 

where t p = t 3 kyr. Meanwhile, ballistic motion of the pulsar away from the SNR's centre gives: 

R p = 1.02V3 pc, (3) 

while we can rewrite Equation (1) to show that: 

r w = 0.01077 1/2 n~ 1/2 v p 1 pc. (4) 

Equations (2), (3) and (4) can be solved simultaneously to give v p = (0.17 - 0.25)^ 1 1 " 04 7j 1 - 21 d§- 3 , n\ = (1.3 - 
2.6)£ 2 1 08 7r 1 - 42 c^ 6 - 6 and t 3 = (93 - lTO)^ 04 ^ 1 - 21 ^ 2 3 . Self-consistency with the requirement that m > n then 
implies d 5 < 1.3i? 2 i 08 ?7~ 2,42 , which combined with 8 d$ > 0.9 (from Hi absorption) and d 5 « 0.9 (from the pulsar's 
dispersion measure) justifies the choices tfe = 1 and E 5 i = 1. 

Thus two different approaches — using either the limits on proper motion, or the morphology of the system — 
both imply a much older and more slowly moving pulsar than results 6 when one assumes t p = r. We note that a 
similar age and velocity to those we have just estimated were derived before the characteristic age of the pulsar was 
known, by assuming n\ = 1 and then solving together the equivalents of Equations (2) and (3) (ref. 7). 

For a pulsar of constant magnetic dipolc moment and moment of inertia, it can be shown that 1 : 

n-l 



t p 



n-l 



1-1^ 



(5) 



where v is the pulsar's initial rotational frequency and n is the "braking index", defined by v — —Kv n (where K is a 
constant), and expected to be n = 3 for a dipole rotating in vacuo. Our lower limit on t p /r from proper motion alone 
implies that n < 1.83, while the value derived from the morphology of the system argues that n < 1.33. These values 
are significantly less than those derived for young pulsars, 18,19 but are similar to n = 1.4 ± 0.2 observed for the Vela 
pulsar, 20 a pulsar which also resembles PSR B1757-24 in its value of r and its tendency to glitch. 21 One alternative 
possibility is that the magnetic field of PSR B1757-24 is growing as a function of time. 22 ' 23 In this case Equation (5) 
does not hold and K is no longer a constant; r is then an approximate time-scale for field growth. 

Further data can support the conclusion that PSR B1757-24 is much older than r = vj2v: multi-epoch VLBI 
observations are expected to show westward motion of the pulsar at the level of ~10 milliarcsec yr _1 , X-ray observa- 
tions of SNR G5.4-1.2 should show emission characteristic of an old rather than young SNR, and a braking index for 
PSR B1757-24 measured from timing data would be expected to lie in the range predicted above. The low braking 
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index of the Vela pulsar, 20 discrepant pulsar/SNR ages in other pulsar/SNR assocations, 23 ' 24 and the fact that many 
pulsars with small characteristic ages have no associated SNR, 2 ' 3 all suggest that PSR B1757-24 is not a statistical 
anomaly. If other pulsars are indeed older than they seem, our understanding of pulsar velocities, asymmetries in 
supernova explosions, the fraction of supernovae that produce pulsars and the physics of neutron star structure and 
cooling must all be reconsidered. 12 ' 25 ~ 27 
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Table 1. Very Large Array observations of the pulsar-powered nebula G5. 27-0. 90 used for proper motion measure- 
ments. Both observations used the same pointing centre, namely (J2000) Right Ascension (RA) 18 h 01 m 00.035 s , 
Declination (Dec) -24°51'27.260". Absolute flux densities were determined using observations of the source 3C 286, 
while antenna gains were calibrated using observations every 15-20 min of the extragalactic source PMN J1751-2524, 
displaced 2?1 from G5. 27-0. 90. The resulting images were formed using natural weighting and multi-frequency syn- 
thesis with 100 mas pixels, deconvolved using a maximum entropy algorithm, and then smoothed with a circular 
Gaussian of FWHM Of! 90. 



Date Observed 
Array Configuration 
On-source integration time (h) 
Centre Frequency (GHz) 
Bandwidth (MHz) 



Epoch 1 Epoch 2 

1993 Feb. 02 1999 Oct. 23 



BnA BnA 

3.3 2.8 

8.44 8.46 

75 100 
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RA (J2000) 



Figure 1. Radio emission from the supernova remnant G5.4-1.2. The image corresponds to archival VLA data at 
an observing wavelength of 90 cm; contours are at levels of 10, 25, 50, 100, 150 and 200 mJy beam -1 , and the peak 
intensity is 150 mJy beam -1 . The resolution of the image is 60" x 45". The pulsar (whose position is marked) is at 
the head of the compact nebula G5. 27-0. 90 to the west of the supernova remnant. 
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Figure 2. Radio emission from the western tip of the pulsar-powered nebula G5. 27-0. 90. The upper panel represents 
the 1999 image of G5. 27-0. 90, with corresponding contour levels at 0.75, 1.5 and 2.25 mJy beam -1 . The ellipse near 
the head of the nebula marks the position of PSR B1757-24, (J2000) RA 18 h 01 m 00.023(2) s , Dec -24°51'27.15(5)", 
which was determined by 20 cm A-array observations in Aug 1999 (uncertainties in the position of the pulsar are 
a factor of five smaller than the extent of the ellipse). The lower panel shows the difference between the 1999 and 
1993 images, on which are superimposed the same contour levels as in the upper panel. To compare the two epochs 
in a quantitative fashion, we took the deconvolved model from each epoch, computed its two-dimensional Fourier 
transform, and then sampled this transform with the transfer function corresponding to the u — v coverage of the 
1993 data. Assuming fidelity in the original deconvolution, this results in two data-sets which have identical u — v 
coverage and sensitivity, and differ only in the actual brightness distribution between the two epochs. We searched 
for proper motion between these two images using the MIRIAD task IMDIFF, 28 which finds the shift in x and y which 
minimises 29 the rms noise of the difference map between the two data-sets. Shifts corresponding to a non-integral 
number of pixels were computed using cubic convolution interpolation. 30 To estimate the uncertainties associated 
with this approach, we repeated this comparison 40 times, each time adding Gaussian noise of rms 30 /xJy beam -1 
(corresponding to the sensitivity of the Epoch 1 data) to one of the images. 



